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During  this  period  the  following  scientists  and  students  were  supported, 
in  addition  to  the  principal  investigator: 


G.  Beylkin 


Post-Doctoral  Scientist 


K.  Golden 


Graduate  Student 


B.  LeMesurier  Graduate  Student 


The  following  papers  have  appeared,  have  been  submitted  for  publication, 
or  are  in  the  final  stages  of  preparation: 


K.  Golden  and  G.  Papanicolaou.  Bounds  for  effective  parameters  of 

heterogeneous  media  by  analytic  continuation.  Comm.  Math.  Phys.  j)0  (1983), 
pp.  473-491. 


G.  Beylkin.  Gaussian  beams  and  representation  of  the  Green's  function. 
Submitted  to  Wave  Motion. 


D.  McLaughlin,  G.  Papanicolaou  and  0.  Pironneau.  Convection  of  Microstructure. 
Submitted  to  SIAM  J.  on  Applied  Mathematics. 


G.  Papanicolaou.  Diffusions  and  random  walks  in  random  media.  In 

Mathematics  and  Physics  of  Disordered  Media,  Springer  Lecture  Notes  in 
Mathematics  it  1035,  1983. 


R.  Figari,  E.  Orlandi  and  G.  Papanicolaou.  Diffusive  behavior  of  a  random 
walk  in  a  random  medium.  Proceedings  of  Kyoto  Conference,  1982,  N.  Ikeda, 
editor,  to  appear. 


G.  Papanicolaou.  Macroscopic  properties  of  composites,  bubbly  fluids, 
suspensions  and  related  problems.  Eyrolles,  France,  1984,  to  appear. 

(These  are  the  lecture  notes  for  a  series  of  13  lectures  given  in  July  1983.) 


D.  Dawson  and  G.  Papanicolaou.  Random  wave  processes.  In  preparation. 


Remarks  on  the  research  effort 


Beylkin,  who  is  now  permanently  at  Schlumberger-Doll  Research, 
Ridgefield,  Conn.,  worked  on  the  use  of  Gaussian  beams  to  represent  wave 
fields  at  high  frequencies  in  a  manner  that  is  insensitive  to  singularities 
such  as  caustics.  Several  investigators  have  been  interested  in  this  subject 
recently  (at  Exxon  in  Houston  and  in  New  Jersey,  at  Schlumberger  and  else¬ 
where),  aJ though  there  is  no  ciear  consensus  yet  regarding  the  usefulness 
of  beams  relative  to  more  conventional  geometrical  optics  methods.  One  of 
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the  principal  objectives  of  Beylkin  was  to  find  how  the  beam  representation 
could  help  in  the  analysis  of  random  media.  The  presence  of  singularities 
is  a  very  serious  limiting  factor  there  (recent  work  of  B.  White  at  Exxon 
has  analyzed  wave  fields  at  high  frequency  in  a  random  medium  up  to  the 
onset  of  caustic  formation). 

Papanicolaou  and  Golden  continued  their  work  on  the  analytic  continuation 
method.  Golden  is  writing  his  dissertation  now  and  will  obtain  his  degree 
in  June  1984.  A  complete  description  of  progress  in  this  area  will  be  given 
later.  The  enclosed  paper  gives  an  idea  of  what  had  been  acomplished  until 
last  February. 

B.  LeMesurier  has  been  working  on  the  focusing  singularity  for  the 
nonlinear  Schrodinger  equation,  exclusively  by  numerical  computations  up 
to  now.  The  analysis  of  the  focusing  singularity  is  our  principal  concern 
now.  We  want  to  develop  accurate  and  well  structured  computational  schemes 
for  analyzing  singularities  in  nonlinear  beams.  As  our  analytical  work  has 
shown  (with  Weinstein  and  McLaughlin  last  year;  the  final  write-up  on  this 
is  still  in  preparation),  the  focusing  singularity  requires  deep  and  difficult 
analysis  for  its  solution.  It  is  also  a  basic  problem  that  arises  in  many 
applications  in  nonlinear  wave  propagation. 

Some  work  on  random  media  is  reported  in  the  enclosed  papers.  This  is 
work  that  was  done  earlier  and  is  continuing  at  a  much  reduced  pace  now. 

The  works  with  Dawson  and  with  McLaughlin  and  Pironneau  were  discussed 
earlier.  The  paper  with  McLaughlin  and  Pironneau  is  included.  We  want  to 
do  much  more  numerical  work  in  connection  with  this  problem  and  we  will 
move  in  this  direction  in  the  near  future.  The  paper  with  Dawson  is  in 
preparation.  (A  copy  of  the  preliminary  version  was  sent  earlier.) 
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A  SUMMARY  OF  WORK  THIS  SUMMER  ON  THE 


NON-LINEAR  SCHRODINGER  EQUATION 

B.  LeMesurier 


The  equation  is 

2  2 
iufc  +  Au  +  |u|  u  =  0  on  [0,2tt] 

with  periodic  boundary  conditions. 

The  work  done  was  almost  entirely  numerical,  although  some 
attempts  were  made  to  extend  results  from  the  case  of  u  in 

2 

Schwartz  class  on  J*  . 

Some  previous  work  of  this  nature  has  been  done  by  P.L.  Sulem, 
Cl  Sulen  and  A.  Patera  [1] ,  giving  at  least  one  example  of  initial 
data  that  seems  to  evolve  to  "blow  up,"  that  is  supjuj  -*■  00  in 
finite  time.  I  have  used  their  algorithm,  as  follows 

Writing  u1?  .  -  M  <  j,k  <  M-l  for  the  complex  Fourier 

] ,  k  —  — 

components  of  u  values  on  a  grid  2M  * 2M  at  time  step  n  , 
the  partially  implicit  scheme 


i(un+l-un) 


1 , .  n+1  ,,  n> 
•^-{Au  +  Au  ) 


,  n-1 1 2  n-1 
,u  u 


gives  the  Fourier  transformed  equation 


,~n+l  ~n  ,  ,  •  2 . ,  2  wf-n+l  ,~n  .  ,-,/i  ni  2  n. 

(uj,k'uj,k)  -  (3  +k  >(uj,k  +  uj,k>  +3(iu  1  u  > 


(lun-1l2un-'l 


I  have  worked  with  M  a  power  of  two,  usually  32,  to  allow  greatest 
speed  in  computation  of  Fast  Fourier  Transforms. 

Algorithms  were  written  for  these  as  this  turned  out  to  be 
faster  than  the  available  packaged  code  (N.C.A.R.)  for  general  M. 


Input  Parameters 

In  addition  to  an  array  of  .initial  values  for  u  ,  there  are 

dt  the  time  step  size 

ttime  the  maximum  time  to  which  solution  will 

be  continued 

umax  a  maximum  value  for  max  u 

The  solver  stops  when  either  the  ttime  or  umax 
limit  is  attained. 
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graphics  output  is  generated 
every  IPF  time-steps 


printed  output  as  specified  below 
is  generated  every  ISF  time-steps 


Typically  ISF  <<  IPF 


Output 


The  graphics  output  currently  consists  of  an  array  of  values 
of  ju|  ,  with  boundary  values  repeated  at  both  boundaries.  This 
is  input  to  contour  plotting  and  other  programs. 


The  current  time  is  also  output.  {Other  graphics  output, 
in  particular  phase  information,  is  now  being  added.) 


The  printed  data  consists  of  a  list  of  the  above  input 
parameters,  and  the  following  (Additions  are  frequently  made  to 
this  list  as  new  questions  arise) 


=  ’  |u| 


H  =  j  ( | Vu | 2  -  ||uj4) 


These  are  conserved  by  the  P.D.E.  but  not  the  algorithm,  so 
measure  accuracy  of  the  solver. 


max|u|  "I  These  three  should  all  "blow-up"  for 
1  suitable  initial  data. 


{  |u| 


The  remaining  Quantities  are  "localized";  the  integrals  are  taken 
over  a  disc  in  the  center  of  the  period.  This  part  of  the  output 
was  added  after  summer,  principally  to  study  the  variance  of  part 
of  the  data,  which  in  the  Schwartz  class  data  case  is  the  key  to 
proving  blow-up. 


Thus  a  new  input  parameter  RL  is  added,  and 


2  2  2 

B  =  {  (x,y)  :  jx-ir[  +jy-uj  <  (2nRL)  }  in  the  following  list  of  output 


quantities 


VlccT  J 


r2  ♦ 


*  2 

Various  factors  of  4 nr  etc.  are  omitted, 


c*  **»  *  I •*.  ♦  *,  •*,  »*,  •*  .*  .*  ,*  .•  .*  .•  ,*  *,•  t*  -  •  .«  ♦  • 
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Observations 


Initial  data  has  mostly  been  chosen  with  a  "bump"  shape: 
one  or  two  maxima  near  the  period's  center,  dying  off  to  small 
values  (and  small  gradients)  near  the  boundary.  All  have  been 
real  valued,  and  most  have  had  H  <  0  ,  a  condition  guaranteeing 
blowup  in  the  Schwartz  class  (i.e.  all  space)  case. 

The  results  has  been,  consistently,  growing  and  narrowing 
of  the  bump  at  an  accelerating  rate,  as  expected  from  the  all 
space  case.  However  limitations  on  resolution  only  allow  a 
growth  in  max|u[  by  a  factor  of  about  3  before  computational 
errors  become  substantial:  a  finer  grid  will  be  needed  in  later 
work . 


As  one  hope  is  that  the  evolution  of  a  bump  to  blow-up 
is  mostly  local  (not  affected  much  by  values  away  from  the  bump 
so  long  as  they  are  small  or  flat,  say)  experiments  have  been 
done  where  the  central  bump  is  left  unchanged,  but  the  values 
outside  are  set  to  various  different  constants  (smoothly  connected 
to  the  bump) . 

The  results  are  promising:  the  rate  and  nature  of  the  blow-up 
changes  little,  and  the  near-constant  outer  region  evolves  much 
as  if  the  bump  were  not  present.  That  is,  it 
retains  constancy  in  space  and  the  same  absolute  value,  to  a 
fair  degree.  A  second  localization  test  was  the  use  of  tv/o 
separated  bumps,  as  below. 


The  evolution  was  compared  to  the 
evolution  of  data  with  just  one  or 
the  other  bump. 


Again  the  results  were  pleasing: 
the  two  bumps  seem  to  evolve  in  ignorance 
of  each  other. 


3 


rt  TT  Cvco<2  dnod<rw 


Later  work  (after  the  summer)  has  studied  the  localized 
variance  of  a  bump,  as  mentioned  above.  The  motivation  is  the 
result 


d2V 


=  8H 


V  = 


df 


,  2 


u(x,t)j  |xj  dx  for  the  all-space  case 


Local  variances  have  been  studied  with  the  cutoff  well  down  the 
shoulder  of  a  bump,  with  and  without  radial  symmetry,  and  with 
the  value  outside  the  cutoff  not  necessarily  going  to  zero.  (The 
value  near  the  boundary  has  always  been  a  constant,  initially.) 


The  local  variance  has  consistently  evolved  with 


V 


dt 


loc 


loc 


In  the  case  of  localization  containing  the  support  of  the  initial 
data,  the  ratio  seems  to  be  a  little  less.  This  is  reasonable 
as  u(x,o)  =  C  ,  a  constant  gives  v^oc  =  constant,  H  <  0  . 

As  yet,  V.  has  not  got  near  0  in  the  runs,  and 

it  is  not  clear  whether  or  not  it  would  approach  close  to  0 
before  blow-up  occurs.  This  question  awaits  higher  resolution 
computations . 


II]  P.  L.  Sulem,  C.  Sulem  and  A.  Patera,  "Numerical  Simulation 
of  Singular  Solutions  to  the  Two-dimensional  Cubic 
Schrodinger  Equation."  In  the  thesis  of  C.  Sulem, 
L'Universit4  Paris-Nord. 


Thesis  of  Ken  Golden.  Outline. 


1 .  The  Analytic  Continuation  Method  for  Two-Component  Media 

Let  (fi,F,P)  be  a  probability  space  and  e(x, a>)  be  a  stationary  random  field 
on  x  €  R  and  <o  €  H.  For  two-component  media  let  £(&)  =  +  e2^2^ 

be  the  dielectric  constant  at  the  origin  in  R^  of  the  realization  co  £  £1  of 
the  material,  where  =  1  if  medium  1  occupies  the  origin  and  =  0  other¬ 
wise,  similarly  for  X2  •  The  effective  dielectric  constant  of  the  composite 
may  be  defined  as 


P(d£o)£(tj)E(a3)  *e 


where  E  is  the  electric  field  at  the  origin  with  average  e^  ,  a  unit  vector 
in  the  It1"*1  direction.  It  is  useful  to  consider  the  function  F(s)  =  1  -  m(h) 
where  m(h)  =  and  5  =  l/(l-h)  .  Maxwell's  equations  yield  a 

resolvent  expression  for  E  so  that 

-1-* 

F(s)  =  P(do)x1(s  + rxx)  ek*ek  (2) 

n 

where  F  =  V(-A)  .  The  spectral  theorem  in  an  appropriate  Hilbert  space 

then  gives 

K.)  -  f  ,  (3) 

j  s  —  Z 

0 

where  y  is  the  positive  Borel  measure  on  [0,1)  arising  from  the  family  of  projec¬ 
tions  associated  with  the  self  adjoint  operator  Tx^  •  The  representation  (3) 
can  also  be  proved  by  noting  that  Im(-F)  >  0  when  Im(s)  >  0  and  that  F  is 
analytic  off  [0,1]  with  F(°°)  =  0.  Then  a  general  theorem  in  function  theory 
gives  (3). 

For  1  s  |  >  1  we  can  expand  (2)  about  a  homogeneous  medium  (h  =  1  or  s  *  <*>)  , 


(  rx, 

F(s)  =  P(du) 

J  i.  s 


cl'kkxl  X1 

2 

s 


Vkk;/r  kkxi 
3 


where  F. 


Equating  (4)  to  the 


expansion  of  (3) 


2 


determines  the  moments  of  y  .  Then  (3)  provides  the  analytic  continuation  of 

(4)  from  1 s [  >  1  to  the  full  domain  of  analyticity  of  F(s).  If  (4)  is  assumed 

known  only  to  first  order,  then  only  the  mass  of  y  is  fixed,  where  p^ 

is  the  volume  fraction  of  medium  1.  Under  this  assumption  we  obtain  for  fixed 

s  ^  [0,1)  extremal  values  of  the  linear  functional  Fg(y)  by  evaluating  (3) 

with  Dirac  point  masses,  since  they  are  the  extreme  points  of  the  set  of  positive 

Borel  measures  on  [0,1)  with  fixed  mass.  The  resulting  bounds  restrict  F  (y) 

s 

to  a  region  in  the  complex  plane  bounded  by  two  circular  arcs,  which  give  for 
real  s  the  classical  Wiener  bounds  (p^/e^  +  ^2^2^  ^  —  e*  —  Piei  +  P2e2  ' 
the  material  is  further  assumed  to  be  statistically  isotropic,  then  it  is  possible 
to  calculate  the  second  term  in  (4).  This  F  known  to  second  order  can  be 
transformed  to  a  function  of  the  type  (3)  known  only  to  first  order,  the  extremi- 
zation  of  which  gives  the  Hashin-Shtrikman  bounds,  which  were  first  found  using 
a  variational  principle.  The  transformation  procedure  works  to  any  order  in  (4). 


2.  Extension  to  Multicomponent  Media 


As  Bergman  [6]  and  Milton  [7]  point  out,  it  is  not  clear  from  the 
two-component  case  how  to  extend  the  analytic  method  to  three  or  more  components. 

We  now  describe  such  an  extension.  For  simplicity  we  consider  only  three-component 
materials.  The  effective  parameter  e*  is  then  an  analytic  function  of  two 
complex  variables  h^  =  e^/e^  and  ^2 = e2^e3  with  Im  e*  >  0  on 

{im  h^  >  0}  *  ilm  h2  > 0)  .  The  key  step  in  the  extension  is  to  obtain  a  represen¬ 
tation  formula  analogous  to  (3).  Using  Cauchy's  theorem  in  several  variables  we 
have  proven  the  following  for  the  polydisc  =  {jC^j  <  l}  *  U^i  <  :  for 

f ( ^i , ^2 )  to  b®  analytic  with  non-negative  real  part  in  it  is  necessary 

and  sufficient  that  f  may  be  represented  as 


f(C1,;2)  =  iv(0,0)  +~  [  l  (H1H9  +  H1  +  H9  -  l)dy(t,,t0)  , 


2  J  J  ^1«2  +  H1  +  H2 


'1*  2 


it .  it . 

where  H.  =  (e  J+C.)  /  (e  , 


q  ’ 


(5) 


j  =  1,2,  v  =  Imf,  and  u  G  M 


3 


>1^  =  ^positive  Borel  measures  of  mass  q  on  the  torus 


satisfying  the  condition 


r  r  i(nt1+mt2) 


dy(t^, t9)  =  0 


(6) 


when  nm  <  0 


I 

I  ■ 


We  evaluate  (5)  with  dy(t^,t9)  =  * dt2 + dt^ *  dv9(t2)  ,  where  and 

v2  are  positive  measures  on  the  circle.  Let  U  =  {Ira  s1  > 0}  *  {im  > 0}  , 

where  each  factor  U  is  conformally  equivalent  to  the  disc  D  ,  with  s  =  1  /  (1-h^) 

and  s2  =  1  /  (l-h2>.  Now,  a  function  K(s^,s2)  analytic  on 

U2UU2U  (([0,1)  *1R)  U  (  TR*  [0,l)))c  with  Im  K  >  0  on  U2  and  K(°°,°°)  =  0  , 

2 

where  denotes  complex  conjugation  and  "c"  denotes  complementation  in  ]R  , 
has  the  representation 


K(s1,s2) 


|  dy  (z  )  j  dy  (z9) 

'  — —  +  •  “ 


0 


*r*i  j,  z2-s2 


(7) 


where  y^  and  y2  are  positive  measures  arising  from  and  v2  .  If 

K(s^,s2>  =  -F(s^,s2)  then  the  expansion  of  F(s^,s2)  corresponding  to  (4) 

forces  the  masses  of  y^  and  y2  to  be  p^  and  p2  ,  the  volume  fractions 

of  media  1  and  2,  respectively.  If  y^  and  y2  are  point  masses  at 

z.  =  z9  =  0  then  y  on  the  torus  is  an  extreme  point  of  M  .  Such  y  give 

l  l  *>l+f>2 

the  real  Wiener  bounds.  When  the  material  is  statistically  isotropic  F  has  a 
non-zero  second  order  crossterm  in  its  expansion,  so  that  (7)  is  not  applicable. 

However,  the  function  G  =  F  /  (1  -  ^  F)  has  the  expansion 


Pi  Py  Pi  P?  0 
G(s  ,s9)  =  —  +  —  +  — 9+— 9+  ---  +  ... 
12  s.  s_  ,  2  2  s.s9 

1  2  ds^  ds9  1  2 


point  masses  at  z^  =  z2  =  1/d  gives  the 


.  Letting  K  =  -G  with 

real  Hashin-Shtriknau  bounds 


and 

for 


y2 


three-component  materials. 


Currently  we  are  working  on  proving  for  the  first  time  complex  versions 

of  the  Wiener  and  Hashin  -  Shtrikman  bounds  for  three-component  materials  which 

restrict  F  to  regions  bounded  by  circular  arcs.  In  our  formalism,  establishing 

these  bounds  and  the  previous  real  ones  depends  on  showing  that  for  fixed 
2 

((,.*£,,)  £  D  ,  the  extremal  values  of  the  linear  functional  f,  (y)  on  M. 

1  A  sl’^2 

are  attained  when  dy  =  *  dt£  or  dy  =  dt^  *  6t*  ,  where  t*,t*  e  [0,2tt)  . 

We  have  evidence  that  such  a  statement  is  true  and  are  working  to  prove  it. 

These  matters  will  be  contained  in  the  author's  Ph.D.  thesis. 

Underlying  a  full  understanding  of  the  three-component  problem  is  the 

characterization  of  the  set  of  extreme  points  of  M^.  Rudin  [5]  and  McDonald  [8] 

have  given  examples  of  these  extreme  points  but  a  full  characterization  of 

them  is  still  unknown.  It  is  not  even  known  if  all  the  extreme  points  are 

2 

singular  with  respect  to  the  Haar  measure  on  T  .  We  propose  to  analyze  the 
extremal  structure  of  .  One  way  of  gaining  insight  is  to  consider  the  set 


of  positive  measures  satisfying  the  Fourier  condition  in  (6)  living 

2  2  2  r  0  1  n- 1  ■» 

not  on  T  ,  but  on  a  discrece  model  of  T  ,  namelv  T  =  iu  ,c»i  /  * 

n*n 


w  r  0  1  ri-1-,  i2ir/n 

x  Igj  ,<u  , . . .  1  ,  a)  =  e 


For  a  measure  d„  =  >  a. .5  .  .on 

.  .  n  1J  1  J 

1,J=0  J  OJ  ,0JJ 


2  2  1 
T  ,  condition  (6)  gives  (n-1)  +1  equations  in  the  n“  unknowns  a..  . 
nxn  J-J 

We  have  an  algorithm  which  enumerates  all  the  extreme  points  of  n  ,  which 

2  2  1 
,  _  ,  _,n  -(n-1)  -1  ,,  -  ,  ,  ,  .  ,  ,,nxn 

is  a  convex  compact  subset  of  3R  .  Hopetully,  analysis  of 

will  shed  some  light  on  the  structure  of  .  We  should  also  mention  that 

McDonald  [8]  has  established  an  isomorphism  between  M  and  a  certain  family 

Q  of  linear  operators  on  the  disc  algebra.  Analyzing  the  extreme  points  of  Q 

should  give  information  about  those  of  . 

In  [A]  we  rederive  the  classical  Uashin-Shtr ikman  bounds  for  three-component 
materials  using  analytic  continuation.  Milton  [7]  has  shown  that  for  certain 
volume  fraction  regimes  these  bounds  are  not  optimal,  i.e.  there  does  not 
exist  a  material  that  attains  them.  He  conjectures  tighter  "bounds,"  but  as 
yet  they  remain  unproved.  We  hope  to  settle  this  question  using  the  methods 


5 


we  have  developed.  One  approach  is  to  understand  more  fully  the  fractional 
linear  transformations  of  F(s^,s2)  that  we  use  to  derive  bounds  on  £*  ,  in 
the  hope  of  finding  one  that  might  yield  better  bounds  for  certain  volume 
fractions.  Since  certain  transformations  of  F  induce  a  nonlinear  mapping 
from  M  into  M  ,  a  related  approach  is  to  look  for  extreme  points 

q  q 

different  from  the  ones  used. in  the  beginning  of  this  section  that  would 
perhaps  give  better  bounds  for  certain  volume  fractions. 

Finally,  in  section  1  we  alluded  to  a  procedure  which  gives  bounds 
incorporating  nC^  order  information  in  (4).  We  hope  to  do  the  same  for  multi- 
component  media  by  using  fractional  linear  transformations. 


3.  Application  to  the  Conductivity  Problem  in  Bond  Percolation 

If  the  heterogeneous  "material"  is  a  square  lattice  composed  of  bonds  of 
conductivity  with  probability  p^  and  a 2  with  probability  p2  =  1  —  p ^  , 

then  at  percolation  (h=a2/a2  =  0  or  s  =  D  the  expansion  of  (3)  gives 

F(l)  =  U(0)(Pl)  +  h(1)(p])  +  h(2)(px)  +  ...  ,  (8) 

where  F  =  1  -a*lc  ,  a*  is  the  effective  conductivity  of  the  lattice  and  the 
/n  \ 

p  are  the  moments  of  p  .  Expansion  (4)  gives 

p(n)  =  (-l)n  j  P(dco)X1((rX1)nek)-ek  .  (9) 


Bergman  and  Kantor  [9]  realized  that  since  the  bonds  are  independent,  the 
_ _ _ _  , ,(n)  _ i _ _ i  n •  „  _  _  , c _ i _ _ _ ,■  i  „  -  „  ..  . . *  . u  _  . 


moments  p 


can  be  calculated  in  terms  of  polynomials  in  p^  with  coefficients 


composed  of  infinite  series  involving  the  lattice  Green's  function  g(x,y)  with 
x,y  S  .  Independently,  ve  computed  the  first  four  moments  of  p  with 

;j(0)  =  Px  ,  U(1)  =  PlP2/d  and  p(2)  =  ^  'L  (DkL)"g(0,0))2^  A-^p^  -  (DkD~g(0,G;  >pr\ 


l:Mi 


where  D  D.  denotes  the  k, i  second  difference,  and 
k  j 

N 

*  /  \  ,  *•  ,N,  ,  N-k+1 

AMtpl>  <"1>1>  +  J,  Pi 

k“l 


t 


:*•  „ '  ,  • . , .  > .  *  _  ■ 


<m  <' 


<  .  »:  ■  .  tr,*-; 


.  n  ■.  w.  w;  a  v.  vvv 


The  fourth  moment  involves  one  of  the  above-mentioned  infinite  series. 

Note  that  for  d  =  2  the  Keller  relation  m(h)m(l/h)  =  1  determines  the  odd 
moments  of  y  in  terms  of  lower  order  even  moments.  While  we  understand 
the  general  structure  of  the  nt‘1  moment  and  could  with  great  patience  compute, 
say,  the  20c^,  we  have  as  yet  been  unable  to  write  down  a  general  formula 
for  y^n\  To  understand  the  percolation  threshold  one  must  understand  the 
properties  cf  y^  for  large  n.  The  hope  is  that  perhaps  for  d=2  we 
will  be  able  to  extract  enough  information  about  the  y^n^  to  analyze  the 
critical  exponent  and  critical  probability.  We  would  like  to  apply  graph 
theoretic  methods  to  the  problem  as  the  expansion  (A)  is  similar  to  those 
that  arise  in  quantum  field  theory  and  statistical  mechanics.  Another  approach 
is  to  use  the  procedure  alluded  to  at  the  end  of  section  1  to  obtain  tighter 
bounds  on  F(l)  as  more  of  the  moments  of  y  are  known. 
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